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Abstract 

Changes (returns) in stock index prices and exchange rates for currencies are argued, 
based on empirical data, to obey a stable distribution with characteristic exponent 
a < 2 for short sampling intervals and a Gaussian distribution for long sampling 
intervals. In order to explain this phenomenon, an Ehrenfest model with large jumps 
(ELJ) is introduced to explain the empirical density function of price changes for 
both short and long sampling intervals. 
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1 Introduction 



This paper introduces a simple stochastic model, an Ehrenfest model with 
large jumps (EL J), to explain the observed fact that the density function 
of price change fits well to a symmetrical Levy stable distribution near the 
central part of the density function, while it has a truncated tail far from the 
center [7]. The proposed model also explains that it is close to a Gaussian 
distribution for long sampling intervals. 

Changes in stock prices and exchange rates for currencies are one of the best 
examples to which the random walk concept is applicable [5] . We have seen a 
lot of work of related topics since the work of L. Bachelier in 1900. In financial 
theory such as the argument on the Black-Scholes model, the distribution 
of (logarithmic) stock returns is assumed to be Gaussian. It seems that this 
assumption is well accepted for a long sampling interval of empirical data, 
namely those for sampling intervals of several days, a week or longer. 

In the case of a short sampling interval of empirical data, i.e. those for sam- 
pling intervals of several seconds, minutes or hours, the price returns are not 
distributed according to a Gaussian law [2]. The value of a for the empirical 
density function of the S&P500 index estimated for short sampling intervals 
indicates that the density function follows a stable law (a < 2) [7,8]. 

The price changes are therefore said to be distributed according to a stable 
law with a characteristic exponent a < 2 for short sampling intervals, and a 
Gaussian distribution (a = 2) for long sampling intervals [3]. 
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2 Ehrenfest model with large jumps 

The Ehrenfest model [1,4] envisages two boxes + and — , with 2R particles 
distributed in these boxes. A particle is chosen at random and moved from 
one box to the other and the same procedure is repeated. 

We consider a generalized Ehrenfest model in which a V steps of the Ehrenfest 
model take place with probability C/Qi at each step, where a, b > 0, Q > 
1, j = 0, 1, 2, ... and C is a normalization constant defined by: 

For simplicity, a is fixed to be 1 in the following discussion. This generalized 
model becomes identical to the original Ehrenfest model in the limit Q — > oo. 

Suppose that initially there are R + n (-R < n < R) particles in box +, after 
repeating the above procedure s times, there are R + m particles in box +. 
The probability of this event calculated in case of the original Ehrenfest model 
[6] can be described as: 



With the duration of each step for the original Ehrenfest model defined as r, 
consider a particle moving along the x-axis in such a way that at the time s r 
it is located at Am. In the diffusion limit: 




(1) 



where 



(l-z) R - j (l + z) R+j = Y,Ci 3) z k . 



k=0 
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A 2 1 

A^O, r^O, — = D, — >j, sr = t, mA^x, 

2r Rt 



we have 



X2 

lim ^2 P{n\m;s) = J P(x Q \x;t) dt, 



xi <mA<i 



where 



P(x \x;t) 
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We denote the particle number in the box + after s steps as M(s), and define 
the changes as: 

Zas(s) = M(s) - M(s- As), 

where As denotes the step interval. Using the above Z^ s (s), we define the 
empirical density Pas(z) a s ; 



E p As{z)= E 

Zl<Z<Z 2 Zl<Z<Z2 



s last 



(3) 



where s — 1, 2, 3, si ast and 5^ denotes the Kronecker delta. 



Using Eq. (1), the probability density function of the changes Za s (s) for the 
proposed model Pf^ J (z) is expressed as: 



z) 



E E 

m,n i 1 ,i 2 ,...,i As ni,n 2 ,..., 



£ ^P(n|n i; ^)^(ni|n 2 ;6 i2 )x 



C 

x^— P(n As \m; b tAs )5 m _ n , k , 



(4) 



where = 0,1,2,.... Using the diffusion approximation, the above equation 
can be approximated by Eq. (6). 
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~ 0.4 



■f(x)=(1-Bxp(-2x))/(1+exp(-x)) 



t 

Fig. 1. Graph of function f(t) = (1 — e - '). 
Assuming that the initial density function is given by a Gaussian density func- 
tion g(x) = P(0\x; oo), the probability density function of the change Z^tit) 
in the original Ehrenfest model, obtained from the diffusion approximation 
(2), is given by: 
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(5) 



where 5(x) is the Dirac delta function. The variance of the change Z&t(t) is 
given by: 



9 D 

Var(Z*(f)) = — (l-e^"). 
7 



Figure 1 shows the behavior of the function: 
/(*) = 1 - e-«. 
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From this figure, we find that when both At\ and Ati are sufficiently large, 
the difference Var(Z^ tl (t)) — Var(Z^ t2 (t)) remains small, even if Ati and Ati 
are very different. 

The time series of the EL J is interpreted as a time series chosen randomly from 
that of the original Ehrenfest model. Using Eq. (5), the probability density 
function of the ELJ is approximated as: 



pELJ/N \p ° pEhr. / A 

r rAs\ z ) - Qi 1+ ... + i As r T(bn+- + VAs)\ z ) 
il,...,M s ^ 



Qh+-+iA S y 4tt_d(i _ exp{-7 r(b h H h 6 izis )}) 

6XP (~4D1 -exp{- 7 r(^ + + ' ^ 

and the variance as: 



? D C As 

— E 75i7^(l-exp{- 7 r(^ + ••■ + &-)}). (7) 



When 2; is large, the density of ELJ in Eq. (6) will be approximately given by: 



with 



c"^ E ()J jA . (- l! - 

.71,— J As ^ 

where jj of the summation is an integer satisfying the condition 7 r(6 J1 + • • • + 
fr 748 ) S> 1. From the above equation, we find that the density of the proposed 
model decays in proportion to exp(— z 2 ) toward the end of the tail as seen 
in Fig. 3 (a), whereas in Fig. 3 (b), the tail decays in proportion to z~ a ~ x to 
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mid-way when As is small. 



The Levy stable symmetrical density function is given by: 



L(z, As I a, 7) = — / exp(— 7 As q a ) cos(g z) dq, 

TV J 
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for the characteristic exponent a ( < a < 2 ) and positive scale factor 7. 
This reduces to a Gaussian density function when a = 2. The Levy stable 
density is characterized by a fat tail and kurtosis, is consistent with empirical 
data, and satisfies a scaling law. It is one of the most natural generalization 
of the Gaussian distribution which is stable with respect to convolution. 

From Eq. (6), we find the reason that the stable density function a < 2 is 
chosen for the best fitting density function of the proposed model for shorter 
As by the maximum likelihood method that we will explain later. The large 
coefficient of z 2 in the exponential function gives a sharper peak at the origin, 
in contrast, the small coefficient of z 2 affords a flatter peak, broadening the 
density function. Therefore, the density function given by the sum of these two 
types of terms with dissimilar coefficients has a fatter tail than the Gaussian 
density function. 

When As is sufficiently large, the density function on proposed model becomes 
the sum of the Gaussian density functions of almost the same shape. Therefore, 
the density function in the proposed model is given by the Gaussian density 
function for longer sampling intervals. 

The proposed model well describes the changes in the stock indices such as the 
S&P500, TOPIX and currency exchange rates. Let each particle in the box + 
represent a buy stance dealer and each particle in the box — represent a sell 



7 



stance dealer. These dealers will get a lot of information from a newspaper, 
television, radio, business reports etc. on business conditions, weather reports, 
governmental decisions etc. and their view will be changed from buy to sell or 
from sell to buy on the basis of such information. These changes will generate 
a change in the stock price or the exchange rate. We assume that a change 
in the dealer's mind directly causes a change in the price so that M(s) is 
interpreted as the price. 

As the information is spatially and temporally distributed, and any piece of in- 
formation may have a different impact on the dealers, the price will sometimes 
change drastically within minutes, or sometimes only slightly over hours. The 
difference between the progress of time in the ELJ and in the original Ehren- 
fest model describes this phenomenon. The parameters b and Q might thus 
represent the dealers' (price change) sensitivity and frequency of information. 

3 Maximum likelihood method 

The maximum likelihood method [9] is employed to estimate parameters of 
the model. It is assume that a joint density function f(x 1 , x n \9) of the 
random variable (Xi, X n ) is given, where the parameter 9 specifies the 
density function. If the random variables are mutually independent, then the 
joint probability density function of (X 1: ...,X n ) is given as the product of the 
individual density functions of Xi, {i = 1, 2, n). Thus we have 

f(x 1 ,...,x n \e) = f(x 1 \e)f(x 2 \e).--f(x n \e). 

Taking the logarithm on the right-hand side yields the log likelihood: 
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S 

Fig. 2. Trajectory of the number of particles in box + (M(s)) for N = 10000, b = 3, 
Q = 2. 

n 

i(9) = ^ogf(x l \e). 

i=l 

For the present case, the log likelihood is expressed as: 
K a n) = J2 lo & L ( Z As(s), As | a, 7). 



The maximum likelihood estimate of the parameters of the model is obtained 
by maximizing this quantity. 



4 Numerical simulation 



The results of a numerical simulation of the trajectory of M(s) is shown in 
Fig. 2, where N = 10000, 6 = 3, and Q = 2. Big jumps that do not exist in 
the original Ehrenfest model are observed in the proposed model. Because of 
these big jumps, the empirical density Pa s {z) of the new model is a better fit 
to a stable distribution (a < 2), whereas the empirical density Pas{ z ) is given 
by a Gaussian distribution in the original Ehrenfest model. 
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(a) 



(b) 
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Fig. 3. Empirical densities Pas(z) for As = 8, 64, 512 (dotted line). N = 10000, and 
6 = 3, Q = 2. The Lvy stable density estimated using the maximum log likelihood 
method is also shown (solid line), a = 1.310, 1.763, 2.0 and jAs = 9.223, 335.0, 
2418. (a) Semi-log plot, (b) log-log plot. 

In Fig. 3, the empirical density Pas( z ) is shown for different time intervals 
As = 16, 64, 512. These results were obtained for a sequence of about one 
million steps of the ELJ. The density obtained by simulation using the pro- 
posed model is well described by a symmetrical Levy stable density function 
(8). 
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Table 1 



a and ^yAs estimated by the maximum likelihood method for a = 1. 





As 


16 


32 


64 


128 


256 


512 


1024 


N = 


= 10000 




a 


1.386 


1.529 


1.763 


1.958 


2.0 


2.0 


2.0 


6 = 


3, Q = 


2 


As 


22.81 


72.81 


335.0 


1285 


2138 


2418 


2528 


TV = 


= 1000 




a 


1.758 


1.956 


2.0 


2.0 


2.0 


2.0 


2.0 


6 = 


3, Q = 


2 


7 As 


45.69 


139.2 


219.8 


247.1 


249.7 


250.3 


350.2 


iV = 


= 10000 




a 


1.719 


1.756 


1.804 


1.859 


1.932 


1.988 


2.0 


6 = 


2, Q = 


2 


j As 


18.05 


39.28 


90.47 


213.6 


524.9 


1155 


1825 


TV = 


= 10000 




a 


1.713 


1.748 


1.796 


1.871 


1.949 


2.0 


2.0 


6 = 


3, Q = 


3 


*yAs 


15.80 


34.18 


7813 


196.6 


493.7 


1080 


1685 



The parameters a and 'jAs of Eq. (8) estimated using the maximum likelihood 
method for different time intervals are shown in Table 1. In the model, a 
increase to 2.0 with increasing As for different values of N, Q and b, as can 
be clearly understood from Table 1. Here we used the results of numerical 
simulation and considered the numerical values for a and 'yAt for several 
values on N, b or Q in Table 1. However, the probability density function of 
the proposed model can be approximately given by Eq. (6). 

Consider a density function given by a combination of an infinite number of 
Gaussian density functions: 
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It is easy to confirm from the knowledge of domains of attraction [4] that the 
density function of the above equation belongs to the domain of attraction of 
the symmetrical Levy stable density with characteristic exponent given by: 



The density functions given by Eqs. (6) and (9) are identical under a first order 
approximation in the limit 7^0 and As = 1. Therefore, the symmetrical 
Levy stable density with characteristic exponent given by Eq. (10) fits the 
density function (6) well as far as the central part of the density function of 
the ELJ is concerned. 

5 Exchange rate 

In ordered to discuss the empirical density function, we should redefine the 
empirical density function as a continuance version. Price changes (returns) 
are defined by: 



where S(t) is the price at time t and At is a time interval. The empirical 
density function P A t{ z ) is then defined by: 



where 8(z) is the Dirac delta function and t takes nAt', n — 1, 2,...,ni ast . 

The empirical density of the exchange rate between the US dollar and the 
Japanese yen over one year (Jan. '94-Dec. '94) [10] is shown in Fig. 4, along 



Z At (t)=S(t)-S(t-At), 




J25(z-Z At (t)) 
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Fig. 4. Empirical density functions PAt(z) for sampling intervals of At = 2 15 , 2 18 , 
2 21 seconds (dots) and the corresponding stable density functions (lines, a = 1.451, 
1.588, 1.900 and jAt = 23.86, 184.5, 4276) estimated by the maximum log likelihood. 

Table 2 

a and jAt estimated by the maximum log likelihood for the empirical density of 
the exchange rate [10]. 



At 


2 14 


2 15 


2 16 


2 17 


2 18 


2 19 


2 20 


2 21 


a 


1.436 


1.451 


1.441 


1.483 


1.588 


1.666 


1.840 


1.900 


-yAt 


13.84 


24.86 


38.95 


75.46 


184.5 


431.75 


1743 


4276 



with the estimated stable density functions. This figure shows that the empir- 
ical density is close to the stable distribution with a characteristic exponent 
a < 2 for short sampling intervals and a Gaussian distribution for long sam- 
pling intervals. Table 2 shows the maximum log likelihood estimate of a and 
'yAt. When the sampling interval is small, the return is distributed according 
to a stable law with characteristic exponent a < 2. In addition, a becomes 
larger with expansion of the sampling interval and approaches 2. Similar re- 
sults were obtained for German stocks [3]. 
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6 Concluding remarks 

This paper introduced a simple stochastic model (ELJ) to understand that 
the observed density function of the price change is approximately given by 
a truncated symmetrical Levy stable distribution for short sampling intervals 
and a Gaussian distribution for long sampling intervals. 

The elastic potential of the Ehrenfest model makes the tail of the density func- 
tion decay rapider than that of the Levy stable density. The empirical density 
function of the price change should not have an infinitely long heavy tail. 
Thus, the empirical density function given by Eq. (6) seems to be reasonable. 
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